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Abstract 

AdS plane waves describe simple backgrounds which are dual to anisotropically excited 
systems with energy fluxes. Upon dimensional reduction, they reduce to hyperscaling violating 
spacetimes: in particular, the AdS^ plane wave is known to exhibit logarithmic behavior of 
the entanglement entropy. In this paper, we carry out an extensive study of the holographic 
entanglement entropy for strip-shaped subsystems in AdS plane wave backgrounds. We find 
that the results depend crucially on whether the strip is parallel or orthogonal to the energy 
current. In the latter case, we show that there is a phenomenon analogous to a phase transition. 



1 Introduction 



The entanglement entropy is a very useful quantity which characterizes the ground state of a 
quantum many-body system. It can be used as a quantum order parameter to distinguish various 
quantum phases. Recently, entanglement entropy has been actively explored both from the string 
theory side and the condensed matter theory side. For example, the properties of entanglement 
entropy of ground states have been studied from the field theoretic approach [J and the holographic 
approach [2]. However, our understanding of quantum entanglement of excited states in quantum 
many-body systems is still limited at present. 

It is well-known that the entanglement entropy correctly measures the amount of quantum 
entanglement only for pure states. On the other hand, if we compute the entanglement entropy 
for a finite temperature system as an example of mixed state, it includes a contribution from 
the thermal entropy in addition to that from the quantum entanglement. Therefore we need 
to look at pure excited states in order to estimate the quantum entanglement by employing the 
entanglement entropy. One such example is an excited state created by a quantum quench, which 
is typically induced by a sudden shift of parameters in a given quantum system. The entanglement 
entropy increases under time evolution after the quantum quench. In two dimensional conformal 
field theories (CFTs), some analytical results have been known [3j. Holographic calculations have 
been done for any dimension [HI!]. In these examples, the systems are excited homogeneously 
and isotropically. Since the holographic description of quantum quenches requires complicated 
numerical calculations of black hole formation in general, it is quite difficult to get analytical 
results on the behavior of entanglement entropjQ. 

The purpose of this paper is to study entanglement entropy in a more tractable setup. Specif- 
ically we consider a CFT with a constant energy flux T^^. This offers us a simple model of 
anisotropically excited states with energy flow. It is holographically dual to an AdS space with 
a plane wave [Tj EJ [9] (similar solutions in the form of shock waves in AdS have been studied 
previously e.g. [lOl [HI [l2]). To define the entanglement entropy we need to specify a subsystem 
which we trace out. We choose the subsystem to be a strip with a finite width. Interestingly we 
will find that the behavior of the holographic entanglement entropy (HEE) crucially depends on 
the direction of the strip. 

After a light-like compactification, as shown in [TJ [HI [9], these AdS plane wave backgrounds 
become gravity duals with hyperscaling violation: the AdS^ plane wave in particular gives a space- 
time exhibiting logarithmic behavior of entanglement entropy. Such a background can be realized 

^Nevertheless, a universal relation which is analogous to thermodynamics has been found recently in [6] for a small 
size subsystem. 
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in effective Einstein-Maxwell-Dilaton theories [13 H [13 [13 [HI [ISl [13 EOl ED [22] and has been 
expected to be dual to non-Fermi hquids [23 [21] (see also [23[Zl[23[3[2Zl[23[23[33[n]for string 
realizations and further discussion). 

The paper is organized as follows. In section 2, we review AdS plane wave backgrounds. In 
section 3, we present results of holographic calculations of entanglement entropy in the AdS plane 
wave for the case where the strip subsystem is parallel with the energy current. In section 3, we 
study the holographic entanglement entropy when the subsystem is orthogonal to the energy current 
and show that there is a behavior analogous to a phase transition. In section 5, we summarize our 
conclusions. 



2 Reviewing AdS Plane Waves 

The gravity /5-form sector of IIB string theory contains the AdS^ plane wave [7] (see also [3 [9]) 
ds'^ = —[-2dx+dx- +^dxf + dr^] + R^Qr'^{dx+f + R^dnl , i?^ ~ gyM^a'^ , (1) 

^ i=l 

as a solution with no other sources, with Q a parameter of dimension (boundary) energy density, 
and dCl'^ being the metric on 5^ (or other Einstein space). Equivalently, the 5-dim part of the 
metric is a solution to Rmn = —^9mn arising in the effective 5-dim gravity system with negative 
cosmological constant: the deformation is a normalizable one. From the holographic stress 
tensor calculation, one obtains 

Thus we see that this spacetime is dual to an excited state in the boundary AA=4 SYM conformal 
field theory with uniform constant lightcone momentum density turned on. These spacetimes 
correspond to a (chiral) wave on the boundary. Imposing a null energy condition, we have ~ 
Q > 0. 

This spacetime ([1]) can be obtained [8] as a "zero temperature" double-scaling limit of boosted 
black D3-branes, using e.g. [31]: to elaborate, consider black D3-branes {i.e. AdS^ Schwarzschild 
spacetimes) with metric ds'^ = r^[— (1 — rQr^)dt'^ + dx"^ + Yli=i ^^i] + ^2(1*^^4^4) and define t = 
1 x^, = 1 with lightcone coordinates x'^: after boosting by A as — t- X^^x^ , we obtain 



ds^ = 



2 



-2dx+dx- + ^{Xdx+ + X-^x-f + dxf 

1=1 



+ ^2(1 -r4r4) ■ 



Note that we have set the AdS radius R to unity here. Now in the double scaling limit ro — )• 0, A — t- 
00, with Q = fixed, ^ reduces to ([T]). (For tq = 0, this is just AdS^y in lightcone coordinates.) 
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Rewriting this in terms of just Q,rQ, gives 



A 2 



-2dx+dx'- + Qr^{dx+ + -^dx~f + ^dxi 



dr"^ 



i=l 

Prom [31], we see that we now have other energy- momentum components also turned on, 

Turning on a small ro about ([T|), this means r++ is dominant while the other components are 
small. In some sense, this is like a large left-moving chiral wave with ~ with a small 
amount of right-moving stuff turned on. This spacetime interpolates between the usual unboosted 
black D3-brane (A = 1, i.e. Q = ^) and the AdS^ plane wave (A — )• oo, with Q fixed). There 
are two nontrivial scales here, Q and tq: for small tq, we expect that physical observables such 
as entanglement entropy are dominated by the AdS plane wave limit, i.e. by Q, with small tq- 
dependent corrections. 

Let us now consider x+-dimensional reduction of ([TJ as in [7] : the 5-dim part of this metric can 
be rewritten (relabelling x^ = t) as, ds^ = R^{-^ + ^^^^^^ + Qr^ {dx+ - -^f) ■ This gives 
the effective (bulk) 3 + 1-dim Einstein metric and exponents 

dsl = + + = 1, z = 3, (6) 



i=l 

dt 



along with an electric gauge field A = —Qp! and scalar e*^ ~ r. Above, we have compared with the 
hyperscaling violating metric in the form ds'^ = r'^^^^' + ^'^^^a'"'"'^'' )) with boundary spatial 

dimension di. This lies in the hyperscaling violating family "0 = di — 1" exhibiting logarithmic 
behavior of entanglement entropy (as we will review later) and has thus been argued to correspond 
to a gravitational dual description of a theory with hidden Fermi surfaces [231124] (see also [25]). It is 
thus interesting to explore these AdS plane waves further, in particular from the higher dimensional 
point of view. 

More generally, we have the (purely gravitational) AdSd+i deformation, which is the AdSd+i 
plane wave, 

ds^ = — i -2dx+dx- + dxi + dr^ + R^Qr^^^{dx+f , (7) 

the Xi being {d — l)-dim (boundary) spatial coordinates, with Q a parameter of energy density in 
d-dimensions. This is a solution to Rmn = —-^9mn, i-&- to gravity with a negative cosmological 
constant. It is also useful to introduce the time and space coordinate (i,Xrf_i) such that 

x± = '-^^ . (8) 
V2 



This can then be used to obtam the spacetime ([7]) as a double-scahng hmit of the AdSd+i black 
brane in Hghtcone coordinates near a double-scahng limit [5], resulting in a near-extremal AdSd-\-i 
plane wave: most of our discussion with the regulated (finite temperature) case will be for the AdS^ 
plane wave. 

Now, dimensionally reducing on the x+-dimension (and relabelling x~ = t) gives the metric 
and exponents 

,2 i?2(i?2Q)W-2) . df^..\ d-2 „ d-2 

i=l 

For the special case of d — 2 = 2, this 6 value lies in the special family "0 = di — 1" as we have seen 
abovqj. From the lower dimensional point of view (the are numerical constants), the d-dim 



action / ^d+i^^ _g{d+i) (^("^+1) _ 2A) dimensionally reduces as 

' dx+d'^x^f^) (i?W - #Ae-2^/('^-2) - #{d(^f - #e2('^-i)<^/('^"2)F2j, (10) 



where the scalar is gd+i^d+i = e^*^, the (purely electric) gauge field is j4 = —^ and the d-dimensional 
metric undergoes a Weyl transformation as g^fj = e'^'^^^'^~'^') gjfi^^K It is straightforward to check that 
the solution ([9]) is consistent with the equations of motion, with the scalar of the form e^*^ = r'^~^. 

There is a more general family of AdS null deformations [9j which include inhomogenous AdS 
plane waves, of the form ds^ = :pz[—2dx~^dx~+Yli=ldx'^ + dr'^]+g+-^-{dx~^)'^ + dQ,l , with g^+{r,Xi) 
spatially varying. This also includes the case 5++ > const = K , sourced by other matter fields. 



which are asymptotic to z = 2 4-dim Lifshitz spacetimes upon x'^'-dimensional reduction [32]. All 
these solutions, including ([T|), have finite curvature invariants everywhere, which are the same as 
those for AdS^ x S^. Furthermore, they all preserve some supersymmetry. Due to the lightlike 
nature, nonzero contractions involving curvature components vanish (this can be checked explicitly 
for low orders) and thus these spacetimes are likely to be a'-exact, somewhat analogous to plane 
waves and AdS^ x S^. This makes them possibly more interesting as string backgrounds. 

It is worth noting the possibility of tidal forces diverging in the deep interior r — t- 00, even 
though curvature invariants are regular and the same as in AdS^ x 5^: this is a general feature 
of plane wave (and in fact any lightlike) spacetimes. This may reflect the fact that from the 
lower dimensional point of view, we have a hyperscaling violating spacetime, which is conformal 
to Lifshitz, and there are curvature singularities arising from the conformal factor (unlike Lifshitz 
spacetimes which only have diverging tidal forces); see e.g. |33j for a recent discussion. From the 
field theory point of view, this is the question of whether turning on uniform r+_|_ density has 



■^Note that the di in the expression 9 = di — 1 is the boundary spatial dimension, while we are discussing AdSd+i 
plane waves with d the boundary spacetime dimension in the higher dimensional description. 
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certain pathologies. From the bulk point of view, this sort of a singularity if it exists is often 
regarded as mild, possibly regulated by finite temperature effects. In this sense, thinking of the 
homogenous AdS plane wave as a zero temperature "chiral" limit of the boosted black brane is 
useful: physically any singularity will be cloaked by the finite temperature horizon. From the dual 
point of view, we are considering a certain boosted limit of thermal states in the CFT: for finite (if 
large) boost, we expect this is well-behaved. 

2.1 Entanglement Entropy for Light-like Subsystems in AdSs Plane Wave 

The entanglement entropy Sa for the subsystem A is defined by Sa = —trpA^ogpA, where pA is 
defined by tracing out the density matrix p for the total system over the subsystem B, which is 
the complement of A, i.e. pA = tr^/o. In any backgrounds of AdS/CFT, we can holographically 
calculate the entanglement entropy Sa from the area of extremal surface ja which ends on the 
boundary of A. An extremal surface is defined by the one whose area functional gets stationary 
under any infinitesimal deformations with a fixed boundary condition. The area is calculated in 
terms of the Einstein frame. If there are several extremal surfaces, we need to pick up the one with 
the minimal area among them. This is called the (covariant) holographic entanglement entropy 
(HEE) d]: 

Area(7A) 



Sa = Mm^^(zEs. 
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where G^^""^' is the Newton constant in the d + 1 dimensional gravity theory we consider. In static 
backgrounds, since we can restrict to a time slice, ^a is reduced to a minimal area surface on that 
slice and leads to the minimal surface prescription [34]. In many parts of this paper we will omit 
the factor — because we are mainly interested in the dependence on the size of the subsystem 
and thus the overall factor is not important for our purpose. 

We now review the logarithmic violation of the area law for entanglement entropy in the AdS^ 
plane wave, using the HEE (llip and finding the area of a bulk minimal surface bounding a subsystem 
A in the shape of a strip in the xi, 2;2-plane. From the point of view of the lower dimensional theory 
obtained by x^-dimensional reduction, it is natural to consider subsystems A that extend along the 
^"'"-direction completely, and the corresponding bulk minimal surface lying on a constant lightcone 
time x~ slice (which corresponds to a constant time slice since x~ = t). 

Consider a strip region in the xi-direction given by —l<xi< /, extending along the 3:2- 
direction: the minimal surface is parametrized by x = x(r), and its area gives the entanglement 
entropy 

Sa = — TFT / / V9++dx^ / = — / dr , (12) 

2G%^ Jo r Jo J r 2G%^ Je r 



where e is the near-boundary cutoff {i.e. the UV cutoff in the field theory). The minimal surface 

dr I 
dx I 

Vl^ - r2, (13) 



has a turning point r* where 4^|r, = 0. The minimal surface then is given by the half circle: 



and thus we can estimate as 



X 



Sa= ^^^^ Llogi . (14) 



Using G)^' = this gives the logarithmic behavior, as expected from the lower dimensional 

theory. We have effectively taken / ^ e, so that the strip width I is macroscopic relative to the UV 
cutoff e in the field theory. When the strip size shrinks to roughly the cutoff, we have a cross-over 
to the UV behaviour in the field theory: in this case, we expect the entanglement entropy for AdS^ 
in lightcone time slicing which vanishes, as vindicated by ()14p for / ~ e. As Q — t- 0, this surface 
degenerates and becomes null, and the corresponding area vanishes. 

Note that this calculation ()14p above arises from just the 5-dim part of the spacetime, so it also 
applies to AdS^ x plane waves dual to various M=l super Yang-Mills theories. 



3 Entanglement Entropy in AdS Plane Wave: Case A 

The purpose of this paper is to give an extensive study of holographic entanglement entropy in the 
AdS plane wave background ([7]) and its regularized one. Especially we will look at the entanglement 
entropy Sa when the subsystem A is given by a (infinitely extended) strip with a finite width / 
at a constant time. Then the entanglement entropy consists of two terms: one is the area law 
divergence [35] and the other is the finite term \36\ [3l] . We expect that the area law divergence 
remains the same as that for the pure AdS dual to the ground state, while the subleading finite 
term will be modified in our excited backgrounds. 

There are two different choices of such systems depending on whether the finite width direction 
is (a) either of {xi,X2, • • •,Xd_2), called case A or (b) Xd-i, called case B (see Fig(T]). We will study 
the case A below and the case B in the next section because in the latter case, a more careful 
analysis is needed. Notice that in case A, the energy current is parallel with the strip, while it is 
orthogonal in case B. 

3.1 Holographic Analysis in Case A 

In case A, the width direction is orthogonal to X-f. and X-_. We can take it to be along xi. Thus 
we need to assume d > 3. Then the subsystem A is specified by 

< xi < /, {x+,x-) = {ay,-/3y), -co < y, X2, X3, • • •, Xd_2 < oo. (15) 
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Figure 1: The two different choices of the strip subsystem A. The gray arrows represents the energy 
current. The left and right are called case A and case B, respectively. 



The extremal surface 7^ is specified by the function x\ = x(r). We define the maximal value of r 

dr 
dx ' 



to be r^, where turns around, changing the sign of '^^ 



The area functional is given by 

Area = ■ ' drr'^'^-^^ ^ {I + x'2)(2/3 + Qr-^a), (16) 

where denotes the (infinite) volume in the y and {x2^ ■ ■ ■,Xd~2) direction. The infinitesimally 
small parameter e represents the UV cut off. 

The extremal surface is specified by the function xi = x{r). The equation of motion leads to 

drX = (17) 

where A is the integration constant. Note that at r = r*, the denominator of the right hand side 
of (|17p vanishes. Thus we find 

1 = 2 dr , . 18 

Jo ^J2p + Qard - A^'^^'^"^) 

Now the area functional is rewritten as 

Area = 2yJayd-2R / -rr • 19) 
L rd-^ + Qard - A^r^i^-i) 

3.2 Behavior of Entanglement Entropy 

The entanglement entropy Sa in a quantum field theory is defined for a subsystem A which is on a 
constant time slice. This corresponds to the choice /3 = a = 1. In this case, the leading divergence 
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is the standard area law 

Sa ~ Oie-^"^-^^). (20) 

The effects due to non-vanishing flux Q turn out to be finite. When Q and / are both large, we 
find I ~ r^, and Qr"^ ~ A'^r'^^'^ ^' >> 1. In this case, we find that the finite part of Sa behaves liktf 

SA\f^nite ~ ±^d-2\/Q " ^'"1 (21) 



The sign in front of (j2ip is + for d < 4 and — for d > 4. In the case d = 4 we need to replace ()2ip 
with 

Sa\ finite ~ 

V2^/Q- log {l-Q-*y (22) 

In general, when I is large, the above result is smaller than that of the thermal entropy, which 
is extensive and is proportional to while it is larger than the finite contribution in CFTs at zero 
temperature, which scales like ~ 0(1^^'^'"'^^). Moreover, the finite contribution is a monotonically 
increasing function of /. This agrees with our intuitive expectation: the active degrees of freedom 
should increase because the energy flow excites the system. 

We can also consider disconnected extremal surfaces deflned by x{r) =const., which extend 
from r = e to the deep IR limit r = oo. However, they always give larger areas and therefore 
do not contribute. When d < 4, it is clear that the IR contributions to the areas diverge as 

drr^~'^/^ = oo. When d > 4, their contributions are flnite, but numerically we could conflrm 
that they are larger than the connected area (jl9p . This situation is similar to the holographic 
calculations in the pure AdS spaces. 

Two more comments are worth making. First, when c? = 3, so that one is dealing with the AdS 4^ 
plane wave and a 2 + 1 dimensional boundary theory, the flnite part of Sa in eq. ()2ip becomes 

Sa\ finite ~ FiVqV/ (23) 

where Vi is the length of the region in the y direction. We see that the log term in d = 4 
is replaced by a power law enhancement going like Vl. While examples of string constructions 
or constructions in gauged supergravity giving rise to hyperscaling violating geometries are well 
known, to our knowledge, the AdS plane wave is the flrst example of a string construction which 
gives rise to such a power law enhancement. 

Second, note that the AdS plane wave solution represents a valid state in any CFT which 
admits a smooth gravity dual since it arises simply from the AdS Schwarzschild solution in the 
limit of inflnite boost. As such, our result for the entanglement entropy in eq. (|2ip is also quite 



^Here we omitted factors like (j+i) are not interested in tiie overall factor which is fixed once the CFT is 



given. 
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general applying, for example in the d = A case to the AA=4 SYM theory, the Klebanov-Witten 
theory from D3-branes on the conifold, and more general N=l superconformal Yang-Mills theories 
dual to AdS^ x (where is the Sasaki-Einstein base space of the Calabi-Yau cone where the 
D3-branes are stacked), and in the d = 3 case to the M2 brane theory, ABJM theory etc. 

3.3 Light-like Limit 

The limit where A gets light-like is given by /3 = 0. This corresponds to a generalization of 
entanglement entropy and its precise definition in quantum field theory can be given by the replica 
method at least formally (analogous to [3]). 

We can also set a = 1 without losing generality. Then we find Z ~ . Thus it is clear that the 
leading term with respect to e(— t- 0) behaves 



as we reviewed in subsection 2.1. Note that the divergence is always smaller than the space-like 
case /3 / 0. 

3.4 Analysis in Regularized AdS Plane Wave 

The above calculations can be repeated for the regulated or near-extremal AdS^ plane wave (j4]), 
with finite boost A. We obtain 




(24) 



Thus, Sa is finite for d < A. For d 



4, Sa is logarithmically divergent [7] 




(25) 




(26) 



givmg 



drX 




(27) 



and thus 




(28) 



For /3 = 0, we have a lightlike 



surface, giving Sa ~ V2^/Q J 



dr 



1 



with a logarithmic 




leading divergence as before. In addition we have a finite subleading piece V2\/Q^ which is cutoff 
independent. For large size, we expect that the surface wraps part of the horizon. 
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For a spacelike surface, we can take a = /3 = 1 as before, and the leading divergence reflects 
the area law. Taking rg small, we expect the extremal surface to be essentially governed by the 
AdS plane wave background, with corresponding entanglement entropy. 

As a check, note that in the limit where Q = i-^- the familiar black D3-brane, taking 

a = /3 = 1, we recover the familiar expression 5*^ ~ V2 f % ,— = for the black D3- 

brane. 



4 Entanglement Entropy in AdS Plane Wave: Case B 

Next we study the case B where the width direction of the strip A is parallel to x^-i- More 
generally we assume that the strip A is defined by 

Ax+ , Ax+ Ax" _ Ax" , , 

< x^ < , < X < , —00 < Xj < 00. 29) 

2~~2 2~~2 

We will denote the regularized length in each of the Xj directions by L(>> I). The dimension d is 
assumed to be d > 2. 

In this background, we again apply the covariant holographic entanglement entropy (jlip to our 
non-static spacetime ([7|). The extremal surface can be specified by 

x~^=x^{r), x~=x~{r). (30) 

The area functional looks like 

Area = 2R'^~^Vd-2 f' -^^1 - 2{drX+){drX-) + Qrd{drX+f, (31) 

where e is the UV cut off as before. 
The equation of motion leads to 

r''"^ Qr'^{drX~^) — drX~ 



A ^1 - 2{drX+){drX~) + QH(9^X+)2 



AB ^1 - 2(a^x+)(9rx-) + QH(9rX+)2 

where A and B are integration constants and we assume A > and B > 0. 
This leads to 

drX~^ = 



^2{d-l) 



jj + Qr'i- 2B 



^ + Qr^-2B 



(32) 



^ - Qr'^-B 

OrX = — =. (33) 
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The turning point r = r^, of the extremal surface is determined by 



In this way, we obtain the relations 



:^ + Q< -2B = 0. 



dr 



(34) 



2 



+ Qr"! - 2B 



j.2(d~l) 

{Qr'^ - B)dr 







(35) 



j.2(d-l) 

In order to take the subsystem A to be space-like, we can assume Ax"^ > and Ax^ < 0. 



Finally, the area of the extremal surface is computed as 



Area = 2/2'^- / ^ 



AB 



(36) 



4.1 Exact Analysis in d = 2 



In the d = 2 case, we can analytically perform the previous integrations. This corresponds to a 
particular limit of the result discussed in [1]. 
We obtain 



Ax 



1 ^ i + ^VQ 

:l0g- 



Ax" = 



This can be solved into 



The length is given by 



A 



tanh(A/QAx+) 



Q 



(37) 



(38) 



Length = R log — + ii loe 



i?log^ +i?loe 



_ 2w, 
WW 

-A" sinh(^Ax^ 



Q 



Finally, the holographic entanglement entropy is found to be 

-2A- sinh( VQAx+) 



Sa = - log 
b 



(39) 



(40) 



We can confirm that this agrees with the result in CFTs by taking the chiral limit of the analysis 
in m. 
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It may be interesting to define the light-like limit (i.e. the limit the subsystem A gets closer to 
light-like) by 

Ax~Ax+ = -e^, Ax+ = finite > 0. (41) 
In this limit, Sa becomes finite: 

^ = = log . log 2 + f (Ax-)^ + O(Q^). (42) 

4.2 The analysis in > 3 

Let us move on to the higher dimensional case. We concentrate on the case where the subsystem 
^ is on a fixed time slice (t =fixed) and thus we require 

Ax+ = - Ax" = (43) 
^/2 

where I is the width of the subsystem A. 

Assume that the width / is very large. The length I gets infinitely large only when the denomi- 
nator ^^(^rf^i) + Qr"^ — 2B develops a double zero at r = r*. This condition and (|34p can be solved 
at this degenerate point as follows 

4(d-l)^*' 

8id-l)d rt^ (44) 
(3(i-2)2 Q 

However, if we plug this in the integral we find that Ax~ gets positively divergent. This means that 
in the limit Ax"*" — t- +oo, we always find Ax~ — t- oo and therefore the subsystem A gets time-like, 
which is not what we want. Thus this tells us that there exists an upper bound for /, which scales 
like Q'^^^'^ for this connected extremal surface. Note also that we need to take i? > 1 to satisfy 
(|43p. We find that the width is vanishing aX B = 1. 

Actually there is another candidate of the extremal surface. This is the disconnected surface 
simply given by x^ =const. We plotted areas as functions of the width 1/^/2 = Ax~^ = —Ax" in 
Figl21 where we subtracted the disconnected surface area from the connected one. The covariant 
holographic entanglement entropy is defined by choosing the one with the minimal area among 
extremal surfaces as in (|lip . Thus our results in FigjS] show that there is a sort of phase transition 
at the width Ic ~ 2.4 for d = 3 and Ic ~ 1.2 for d = A. For I < Ic the connected surfaces are favored, 
while for / > 1^ the disconnected ones are chosen. This behavior is a bit similar to the behavior in 
the confining backgrounds such as AdS solitons [371 EH] • Notice also that we can easily estimate 
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Naively, the existence of the energy flux T_|_+ oc Q > just excites the system in a similar manner 
to finite temperature systems because the system is excited at the energy scale Q^/*^. Indeed, this 
speculation is true in the d = 2 case as can be seen from (|4U|) . However in higher dimensions d> 3, 
as we found here, the result of holographic entanglement entropy is rather different from that in 
the finite temperature system, which has a positive and extensive (finite) contribution to Sa- Our 
results for d > 3, plotted in FiglJl show that there are positive contributions, but they are not 
extensive at all. The presence of the phase transition is a special feature of case B, which does not 
appear in the case A. 

It might be useful to note that in the d = 2 case, the metric ([7]) is equivalent to that of the 
extremal rotating BTZ black hole via a coordinate transformation. However, this is not true in 
higher dimensions. Therefore only d = 2 has non-zero (thermal) entropy and it is very natural that 
we find the extensive behavior in this case. On the other hand, for d > 3, things are different and 
the systems are far from thermal states. 

4.3 An Interpretation of the different results between Case A and Case B 

It is clear that the difference between the case A and case B comes from the direction of energy 
current relative to the strip direction (see Fig{T]). When the energy fiows along the strip (case 
A), the energy does not leak from the strip. Therefore, the system gets simply excited and the 
entanglement entropy increases monotonically with I (see also [6] for the related argument in static 
systems). This intuitive expectation agrees with our result of HEE in section 3. 

However, in case B, the strip extends in the direction orthogonal to the energy flow and thus 
there is a constant energy exchange with adjacent regions. In this non-trivial setup, our holographic 
analysis shows that there is a maximal width Ic ~ Q^^'^ above which Sa gets constant. 

Since the system we consider have the energy flux ~ Q, each excited mode has the wave 



. screening length at finite 



length of order (r++)^/'^ Q^^'^, which is analogous to the therma 
temperature. Thus for a larger subsystem I > Ic, there is no correlatioro or entanglement between 
the deep inside of A and B except the regions near the boundary of A. This explains why Sa in 
the case B does not depend on / when I > 1^. As usual, the large N limit amplifies this phenomena 
and leads to the sharp phase transition. On the other hand, if we consider the correlation in the 
direction orthogonal to the energy current as in case A, then there is no screening effect and thus 
there is no phase transition. 



* Indeed, in the replica method, the calculation of entanglement entropy can be regarded as the correlation function 
of two defects which produce the deficit angles. In two dimensional CFTs, they are reduced to two point functions 
of twisted vertex operators [39] . We can also generalize this into higher dimensions in principle [341 140) . 



13 



Figure 2: Numerical plots of the regularized areas of extremal surfaces as functions of the width 
/ = \/2Ax+ = —^/2/S.x~ of A. We set Q = R = \. The left and right graph corresponds to the 
results for d = 2> and d = 4, respectively. In each graph, the blue curve corresponds to the connected 
surface in the pure AdS space while the red one to the connected one in the AdS plane wave. The 
area is regularized by subtracting the area of the disconnected surface given by =constant. 

4.4 Regularized AdS plane wave: x^-coordinates 

To understand better the phase transition above, let us regulate the AdS plane wave with the 
non-extremal background @ containing a horizon at rg: for concreteness, we consider regulating 
the AdS^ plane wave in terms of the boosted black D3-brane and Case B, with the strip width 
parallel to X3. Then the entanglement entropy is given by 



Sa - V.j ^^^-i^-2(9.xH-)(9.x-) + Qr4(9.xH- + ^5.x-)' . (45) 
The conserved momenta give 

4 

Qr^{drX~^ + ^drX~) — drX~ 



~A 



- 2{drX+){drX-) + Qr^{drX+ + ^drX-f 



r3 drX+ - ^{drX+ + ^dr-X^) 

~AB " 



„4 



(46) 



2{drX+){drX~) + Qr^{drX+ + ^drX-f 



This gives 



dr 



X 



[I - ry) + Qr4(l + ^B)2 - 2B 



r*r^ _ Br^r ^ 



dx+- - ^ (47) 

^ " "^""^ • /^5^(1 - ry) + Qr4(l + ^5)2 _ 2B 
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Then we have 



Qr^ - 5 + ^''1'"^ 
Ax = / dr 



-(1 - ry) + Qr4(l + ^5)2 - 2B 



Ax+ = I dr ' ■ . (48) 



-(1 - r^r^) + Qr4(l + 1^5)2 _ 2B 



The entanglement entropy finally becomes 



Sa-V,I % . (49) 

'^2^2 _ ^252^4^4 _ 25r6 + (1 + ^BYQi 



For To = and (5 = 0, these expressions reduce to those for AdS^. In the limit ro — )■ with Q fixed, 
these expressions can be seen to reduce to ([35]) ([36|) for the Ad^s plane wave. 



For the unboosted black brane Q = ^ with appearing symmetrically, we can set B = 1: 



2 

this then gives drX~ = , ^ ~l = —drX~^, Sa ~ Vo (I* ^ , „ ^ , giving 

= -^{drX~^ + drX~) = and the expected bulk minimal surface on a constant time slice. The 
denominator in these expressions is a factorized limit of those in (j47p . (|48p. (j49p . 

More generally, let us consider the boosted black brane with A2 = ^ ^ 1: the AdS plane wave 
arises in the extreme limit of infinite boost A — t- oo. Then it can be seen that the denominator in 
(j47|l . (j48]l . (j49]l . can always be factorized as in the unboosted case above if 5 = ^ = A^. At this 
factorization point, we have /S.x~ = fn* dr . ^ 2Q/rp) Ax"*" = C* dr . ^ ^ ^ , 

With 5^ ~ ^2/^^7p5i5z^«)(TZ;p) • This gives Ax" = -A2Ax+, .e. = AAx+ + 

A^^Ax" = . Furthermore, we have y/2\drt = drX~ + X'^dr-x^ oc (—B + ^) = 0, i.e. we have 
a constant-t bulk minimal surface. Here t = -^^ +\_ ^ tj^e unboosted time coordinate: recall 

v2 

that the [x^ -f-)- (t, X3)] -coordinate-transformation maps the boosted AdS black brane ([3]) to an 

unboosted brane. Thus the factorization point i? = ^ is simply the case where drt = for the 

''0 

extremal surface, i.e. the surface is at constant-t. This is of course the usual connected constant- 
time extremal surface in the AdS black brane corresponding to the unboosted observer, rewritten 
in x^-coordinates (with the momentum B tuned as above). From the point of view of the boosted 
brane, this is not the most natural: T = is the more natural time coordinate for the boosted 

case. 

With a view to obtaining large values for Ax^ i.e. large size, we note that the parameters 
A, B can be tuned to a location where there is a double zero of the denominator surface, giving 
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Ax^ — )• oo: this happens when 

V{n) = A^B' - A'B\trt - 2Brl + (1 + ^BfQrl' = , 

4 

V'{r,) = -2A^B\trl - QBrl + 5(1 + ^BfQrl = . (50) 

A rough plot of this surface V{r) suggests that the nature of the surface does not change much 
for smaU rg. In particular, we can expand these expressions in powers of Tq to obtain the leading 
corrections 

B^lQrUlQry.^..., . || + -llM + . . . . (51) 



So it appears that the double zero location shifts a little from the AdS^ plane wave values 
but continues to exist with small ro-dependent corrections: Ax^ are positively divergent. This is 
perhaps not surprising since the scale Q^^^ governing the AdS plane wave is widely separated from 
the horizon scale ro in the near-extremal limit with small tq. 

More generally, eliminating the parameter A in (jSOp . we obtain 

ION 



{2Brt - (1 + :^BYQrl 
and a quadratic equation for B which is solved as 



A^B^ = ^— r4 ' (52) 

1 - r^rf 



„ _ /2Q\ 6-7roM + 3roM - 2^3(1 - roM)v/3^^ 

UoV roM(5-3roM) ' ^ ^ 

Expanding about rgrf ~ 0, we see that this root matches the AdS plane wave values with the 
leading ro-corrections (jSip . while the other root i?+ diverges in the AdS plane wave limit (ro — s- 0) 
and is discarded. This admits a real solution and thus the double zero location exists only if the 
discriminant is positive, i.e. 

rt<-,- (54) 

Note that this location is inside the horizon at r^r^^ = 1. 

One might have worried that the phase transition we observed is an artifact due to the singularity 
at the IR limit r — t- oo. Regulating the AdS plane wave geometry by placing an event horizon in 
the IR, we have calculated the holographic entanglement entropy (HEE) to see if there is any 
problem. This regularized AdS plane wave is obtained by boosting the AdS Schwarzschild black 
brane: the x^ -coordinates appear in some ways inconvenient calculationally from the point of view 
of identifying the connected bulk extremal surface for large subsystem size as appropriate for the 
boosted observer. Instead it turns out to be more convenient to use different coordinates as we will 
describe below: we will later map this to the present discussion. 
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4.5 Regularized AdS Plane Wave: other coordinates 

Consider the HEE in the regularized AdS plane wave for the strip subsystem (width /) in the case 
B. This is equivalent to the HEE in the AdS Schwarzschild black hole 

(55, 

r^(l — r^r^) 

for the strip subsystem defined by the interval {AX'^, AX~). They are related by 

Q.r|!, AX-^-^. (5.) 

where = Below we focus on the d = A case. 

By solving the equations for the extremal surface 7^ as before, in the end we find that the size 
of the strip is given by 

rr* ^3 
Ay = 2 dr , (57) 

Je y/(a262_r6)(l_ ^4^4) +52^6' 

r* -br^ 

At = 2 dr =, (58) 

Je (1 - r4r4)y/(a262 - r6)(l - r^) + b'^r^ 

where r,, is defined by 

{aH'-rt){l-ry,) + b\t = 0. (59) 

The area of ja is written as 

Area = 2 ^ • ^ (60) 

Je v^(o262_r6)(l-r4r4) + 62^6 

In order to keep the inside of the square root i.e. the function 

g{r) = {ah'^ - r^){l - r^r^) + b^r\ (61) 

positive, we need to require 



V3 



which is equivalent to g'{r^) < 0. For simplicity, we will set tq = 1 without losing the generality. If 
we define r = r* to be the turning point, we have g{r^) = and therefore we find 



a^b' = ^ '*\ (> 0). (63) 



1 _^4 



In this regularized geometry, there is no disconnected solution as it cannot penetrate the horizon. 
This guarantees completely well-controlled calculations. Then we need to ask if we can smoothly 
take the limit A — )■ 00, where the horizon is pushed into the deep IR and the AdS plan wave is 
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recovered. To see this, it is crucial to check if we can reahze the extremal surface 7^1 corresponds 
to any arbitrary large values of Ay and At with the space-like condition |Ay| > |At|). 

Thus let us study the limits where Ay and At get divergent. They are given by some combina- 
tions of (i) the limit r* — )■ 1, where 7a approaches to the horizon, and (ii) the condition where g{r) 
develops the double zero. Thus first we define an infinitesimally small parameter hy r] = 1 — > 0. 
The double zero condition is equivalent to the condition g' {r^) = and this fixes h in terms of r* 
as follows: 



h 



^(1 



~ 2V^r] + 0{if). 



Therefore we can take the second infinitesimally small parameter e to be 



-b + 



a/3- 



-6-h2\/67?. 



(64) 



(65) 



So we have two infinitesimally small positive parameters r/ > and e > which control the limits. 
Now we define the near horizon coordinate by 



(5 = r* — r. 

Then we can expand g{r) near the turning point r = r^, as follows 

g{r) = cie5 + C26^ + 0{6^), 
ci = 4^/6 + 0(77) + 0(e), 
C2 = 24 + 0(r/) +0(e). 



(66) 



It is also useful to note 



45 + 47] + 0{rf'). 



(67) 
(68) 



We can estimate ()57p and ()58p by using the above approximation. This leads to 



Ay 



At 



d5 



1 



log 



2b / d6 



1 



C2 + \JC2[C2 + cie) 
(2^/6r? - e) 







4{b + r/)\/cie(^ + 028^ \l"ri[c2ri - c\e) 



arctanh 



\/r/(c2 + cie) 



Notice that since 6 > 0, we need to require 



2\/6?? - e > 0. 



(69) 



(70) 



We can think of many different limits, where both e and r/ go to zero. It is always true that Ay 
gets divergent as 



Ay ~ --^loge. 



(71) 
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Now we take the following limit 

e ~ 7/*^ ^ 0, (72) 

with A; > 1. This leads to the divergent At as follows 

1 , e ^ — 1 , ,^ ^ 
At- -log- ~ ——log e, 73 

2 7] 2k 

where we have employed the expansion: arctanh(l — x) ~ — ^ log(j;/2) + 0{x) in the limit x — )■ 0. 
Thus we find 

m V6 k-i ^^^^ 



\Ay\ 2 k 

This shows that we can construct a connected extremal surface which leads to Ay — t- oo (i.e. 
infinitely large width) for any ratio which satisfies 

»<^<# 

This covers all regions we wanted. In this way, we can conclude that the AdS plane wave limit 
A — )• oo is smooth. 

To compare with our previous discussions in the x^-coordinates, we can map the coefficients 
A,B in ([Ml)) and a,b in (l60|): this gives 

^=^' (76) 
and correspondingly for the parameter a. The double zero location in ()53p then maps to 



1 + f ' 
in agreement with ()59p . ()62p . This then recasts (j48p as 

(1^(1 + B ) _ B ) _ ^2 _ (1 _ ,4^4)) 



(77) 



Ax" 



(l-r^r4)^ J (1 + 6)(1 - r4r4)^rj£^(l - r^) - (1 - 62) + r^r^ 



A,.= /i_4^1ll)=/- (78) 

7 (l-rgr)^ y + ^4^4)^a2|i(;L-r4r4)- (1-62) + ^4^4 

We then have 

A2 1 
Ax- = -^(At + Ay) , Ax+ = -^(At- Ay) . (79) 

A constant t surface has Ax~ = — A2Ax''~, as expected. Scaling to the double zero 6—7-0, and 
to the horizon r^r^: — )• 1, with 6 scaling faster as in (|72p . we see that Ax^ acquire log-divergences 
Ax"*" — )• /cilog(..), Ax^ — )• — /i:2A2 log(..) , corresponding to a spacelike surface as long as the 
constants ki , k2 are not identical. The AdS plane wave calculation can be defined as the limit 
rg — )• 0, r^, — )■ oo, rgr^, — )■ 1 with the scaling above: this is effectively equivalent to the disconnected 
surface mentioned earlier. Thus the overall picture of the extremal surfaces appears consistent. 
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5 Conclusions 



AdS plane waves arise as normalizable null deformations of AdS x S spaces [71 [H] : they are dual 
to excited CFT states with lightcone momentum density turned on. These spacetimes are 
likely a'-exact string backgrounds similar to plane waves and AdS 5 x S^. With T++ ~ Q uniform, 
these spacetimes are spatially homogenous and can be obtained as certain limits of boosted AdS 
Schwarzschild black branes [8] . Thus they are among the simplest anisotropic and excited systems 
in holographic setups. The AdS^ plane wave upon x'^-dimensional reduction gives a background 
lying in the hyperscaling violating family with ^^6 = d — 1", exhibiting logarithmic behavior of 
entanglement entropy holographically. In the higher dimensional description, the corresponding 
extremal surface lies on a constant lightcone time x~ slice and extends along the a; "'"-direction. 

In this paper, we have explored holographic entanglement entropy (HEE) for these AdS^+i 
plane waves for strip-shaped (spacelike) subsystems. There are two different cases: case A and case 
B as explained in Fig(TJ In both case, the leading divergent contribution to HEE is the familiar 
area law term. The finite part behaves differently in the two cases. 

When the strip is parallel to the energy flux (case A) , the finite part of HEE is a monotonically 
increasing function of the width / of the strip. Though it is smaller than the extensive thermal 
entropy, it is always greater than the HEE for the pure AdS dual to the ground state. In the AdS^ 
plane wave case, this finite piece grows logarithmically as 5-^*™*'^ V2-v/Q log(/Q^/^) with size I 
and is perhaps a finite reflection of the log-growth in the null case above. It was argued in \23\ [23] 
that hyperscaling violating metrics with ^^9 = d — 1" are dual to systems with Fermi surfaces. 
Indeed, if we regard Q^^* as the scale of fermi momentum ( energy) kp, then 5"^*™**^ agrees with 
the expected behavior in systems with fermi surfaces [23]. It is therefore interesting to obtain a 
deeper understanding of the origin of the logarithmic behavior in the AdS^ plane wave backgrounds 
from a field theory point of view. In particular, one would like to understand if this behaviour is 
due to a Fermi surface or some alternative mechanism. In the AdS^ plane wave case, the finite part 
grows with / like 5^*™*^ V\\fl^fQ^ so that the logarithmic term in the AdS^ case is replaced by 
the power, \fl. To the best of our knowledge, the AdSi^ plane wave is the first known example of a 
construction arising from string theory which gives rise to such a power law enhancement. 

On the other hand, when the strip is orthogonal to the energy flux (case B), we find a phase 
transition such that for / larger than 1^ ~ Q"^^'^, the HEE becomes constant. Using the regulated 
description in terms of boosted AdS Schwarzschild black branes with horizon at rg, we have seen 
that this behaviour persists for small tq: this is expected since the scale Q is more dominant here. 
For general boost, the connected extremal surface can be identified in a certain scaling limit where 
the surface approaches the horizon (and develops a double zero). 
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One possible intuitive explanation of these different behaviors in case A and case B is as follows. 
Since the system we consider has the energy flux ~ Q, each excited mode has the wave length 
of order (r++)^/'^ ~ Q^^'^ along the energy flux direction. Therefore in case B, where the HEE 
measures the correlation in the flux direction, the HEE gets trivial (i.e. constant) for a large width 
because we can regard Ic as the correlation length. This is the reason why we found 
the phase transition phenomenon in our holographic analysis of case B. On the contrary, in the 
direction orthogonal to the flux, the situation is similar to the ground state and the correlation 
length should be infinite. Therefore we do not have any phase transition in case A. It is very 
interesting to analyze the same question from CFT calculations in e.g. free field theories. 
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